c^ 



OPCT-94-1 
hep-lat/ 9407023 



Low- Temperature Series for the Correlation Length in 
■^■. d = 3 Ising Model 

H. ARISUE K. TABATA 

Osaka Prefectural College of Technology Osaka Institute of Technology 

OO i Saiwai-cho, Neyagawa, Osaka 572, Japan Junior College 

Ohmiya, Asahi-ku, Osaka 535, Japan 

> 

m . July, 1994 

(N 

O 

o 



Abstract 



We extend low-temperature series for the second moment of the correlation 
function in c? = 3 simple-cubic Ising model from u^^ to u^^ using finite-lattice 
{iLiI method, and combining with the series for the susceptibility we obtain the low- 

^ ' temperature series for the second- moment correlation length to v?^. An analysis 

of the obtained series by inhomogeneous differential approximants gives critical 
^ . - exponents 2v' + 7' ^ 2.55 and 2v' ^ 1.27. 

X 



1 Introdunction 



The low-temperature series oi d = 3 Ising model or equivalently strong coupling series of 
d = 3 Z2 lattice gauge theory had long been shorter than high-temperature series. Recently 
they have been extended to higher orders |]1], 0, ^, ^ using finite-lattice method P, §|, 01- In 
this method free energy density in the infinite volume limit, for instance, is given by a linear 
combination of the free energy on appropriate finite-size lattices. The algorithm to give the 
coefficients of the linear combination is so simple. In the standard graphical method, it is 
rather difficult to list up all the diagrams completely that contribute to the relevant order of 
the series. The finite-lattice method avoids this problem involved in the standard graphical 
method and enables us to obtain longer series. The method was applied ind = 3 Ising model to 
get the low-temperature series of the true inverse correlation length (which is equivalent to the 
mass gap in lattice gauge theory) |l]], free energy [||, magnetization and zero- field susceptibility 
0, and surface tension (which is equivalent to the string tension in lattice gauge theory) 0]. 

In this paper we apply the method to calculate the low-temperature expansion series for 
the second moment of the correlation function /i2 in (i = 3 simple-cubic Ising model to order 
M^^, extending the previous result of order u^^ by Tarko and Fisher j^. It gives the low- 
temperature series for the second- moment correlation length squared A2 = ^1^ (H to order u'^^, 
when combined with the low-temperature series of the susceptibility. This is longer by five 
terms than the low-temperature series for true correlation length squared A2 that was derived 
from the true inverse correlation-length given in Ref. |I[]. 

In the next section we present the algorithm to obtain low-temperature expansion series 
for fi2 using finite-lattice method. In section 3 the expansion series for /X2 and A2 is given. The 
low-temperature series for true correlation length squared Ag is also listed for comparison. The 
result of series analysis by inhomogeneous differential approximants is described in section 4. 

2 Algorithm of low-temperature expansion 



The second-moment correlation length squared is defined by 



where H2 is the second moment of the correlation function 

/^2 = Jim - Y.i'^i - '^jfisiSj)^, (2) 

with V the lattice volume and Vi = {xi, yi, Zi) the coordinate of the lattice site i, and /xq is the 
zero-th moment of the correlation function or the susceptibility and d is the dimensionality of 
the lattice. Here in this paper we take the lattice spacing a = 1. 

The algorithm to calculate the low-temperature expansion of the second moment ^2 is the 
following. We consider the partition function 

^(/3,/i,r7,7i,72,73) = I]exp(-7^), (3) 

{S^} 

with Hamiltonian 

H = I3Y^ SiSj + J2ih + liXi + 72i/i + 73^i + rir1)si, (4) 

for a three-dimensional rectangular lattice Aq with a volume V = L^ ^ Ly x Lz- We take 
the fixed boundary condition that all the spins outside Aq are aligned, for instance, to be 
{sj = +1}. Then the second moment ^2 is given by 



2 ( d^ d^ &" d^ 

y^'i) V \dhdrj d-jf d'^l d'j'i 



/^2 = Jim_ — ( ^7:^ - ^71 - ^7^ - ^T M^ ^(/^' ^' V, 71, 72, 73) k=r,=7i=72=73=o- (5) 



Let us consider the set {A} of all three-dimensional rectangular sub-lattices of Aq (A C Aq) 
with the volume l^. x ly x l^ and define if of A as 

/ Q2 Q2 Q2 Q2 \ 

where Z{A) is the patition function for A with the fixed boundary condition that all the spins 
outside A is aligned, and define W^ of A recursively as 

W{A) = H{A) - ^ W{A'). (7) 

A'CA 

Note that H{A) and W{A) depend not on the position but on the size Ix, ly, h of A. We know 

H{A,) = Y: WiA). (8) 

ACAo 



Taking the infinite volume limit we obtain 

= , lim Y^ ^ (L,.-/,. + 1)(L, -/, + !) 

= E W{l,,ly,k). (9) 

We can prove that the Taylor expansion of W{h) with respect to m = exp (—4/5) 
includes the contribution from all the clusters of polymers in standard cluster expansion [Q 
that can be embedded into the rectangular lattice A but cannot be embedded into any of its 
rectangular sub-lattice A'(c A). The series expansion of iy(A) starts from the order of u" 
with n = 2{lx + ly + Iz) — "i- So we should take all the finite-size rectangular lattices that 
satisfy 2(Z^ + /j, + Z^) — 3 < A^ for the summation in eq. @ to obtain the expansion series to 
order u^ . 

In practice for the calculation of ■^^hiZ{lx^ly^lz)\h=r)='^i=-y2=i3='^^ fo^ instance, we have 
only to calculate the partition function Z(/5, /i, 77; 71 = 72 = 73 = 0) to order hrju'^ . 



3 Expansion series 



The low-temperature series of ^2 have been obtained to order -u^^ using rectangular lattices 
with cross-section up to 4 x 5. Bhanot's algorithm of calculating the exact partition function 



for finite-size lattices [10| can be applied to the partition function P), in which the necessary 



memory and CPU time are proportional to A^ x 2'"'^'-y and A^ x 2'"=^'" xl^xlyX l^, respectively. 
The calculation was performed on FACOM-VP2600 at Kyoto University Data Processing 
Center and HITAC-S820 at KEK, both of which have about 0.5 Gbyte of main memory and 
1 — 2 Gbyte of extended storage. Total CPU time necessary was about 3 hours. 

The low-temperature series obtained is listed in table |l|, where the coefficients {m„} are 
defined by 

/^2 = E"^«^"' (10) 



Table 1: The low-temperature expansion coefficients {m„} for the second moment of the 
correlation function /i2, {A„} for the second-moment correlation length squared A2, and {A^} 
for the true correlation length squared A2 in li = 3 simple-cubic Ising model. 
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rrin 


A„ 


K 














1 











2 





1 


1 


3 





-1 


-1 


4 





10 


10 


5 


24 


-14 


-14 


6 


-24 


85 


93 


7 


528 


-169 


-201 


8 


-960 


884 


972 


9 


8496 


-2390 


-2510 


10 


-21312 


10212 


10618 


11 


125904 


-30594 


-31250 


12 


-380016 


116134 


118792 


13 


1813416 


-368934 


-378902 


14 


-6046440 


1337519 


1377207 


15 


25675200 


-4435616 


-4547052 


16 


-90096000 


15764526 


16140138 


17 


358481304 


-53464296 


-54602714 


18 


-1289158128 


187665313 




19 


4943015520 


-643021360 




20 


-17962232976 


2242649168 




21 


67393016880 


-7729951680 




22 


-245697661872 


26894409824 




23 


909676085232 


-93043627527 




24 


-3315864327216 






25 


12172005334848 






26 


-44293518847536 







and {A„} by 

A2 = E^n^"- (11) 

n 

The latter has been derived from the expansion series for /i2 obtained here and the low- 
temperature series of the susceptibility obtained by Guttmann and Enting to order v?^ 
( They obtained the series to order v?^ and we have checked using the finite-lattice method 
that their coefficients are correct to order u^^ ). Among the series coefficients obtained here, 
those for ^2 to order u^^ and those for A2 to order v}'^ coincide to those obtained by Tarko 
and Fisher |^ and we have added new 11 terms. 

We also list the low-temperature series for true correlation length squared A2 defined by 
i 

^'^ ^ 2[cosh(m)-l]' ^^^^ 

where m is the true inverse correlation-length ; 

m = - lim ^ log (C(L)O(O)), (13) 

L— +00 L/ 

where 0{€) = J2{i\zi=e} ^i is the summation of all the spin variables in z = i plane. The 
coefficients {A'„} are defined by 

A'2 = E^>"' (14) 

n 

and are derived from the true inverse correlation-length m given in Ref. '^. We note that, 
although the coefficients {A„} differ from {A^} for n > 6, their ratios {X'n/^n} are within a 
range of 1.02 to 1.03 for 11 < n < 17. 



4 Analysis of the series 



In our preliminary analysis, we estimate the critical exponents of /X2 and A2 by inhomoge- 



neous differential approximants |jTT[, in which the approximants to a function f(u) satisfy 

Qm(u)/'(«) + PMfiu) + R^{u) = 0(n^+^^+^+2). (15) 

The approximants are equivalent to Dlog Fade approximants when Rn{u) = 0. 

We first give the results of the analysis for fi2- We plot in fig. |l] its critical exponent 
2z/' + 7' against the critical point Uc consisting of the data given by all the approximants 
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Figure 1: Unbiased estimates of 2v' + 7' versus Uc for /i2; the data includes the estimates 
from all the approximants with N = (p, 0, 1, ■ ■ ■ , 6 and L + M + N = 19 



with L + M + N = 19 and N = (p, 0, 1, ■■•, 6. We find a liner correlation between the 
estimates for the critical point and the exponent. Linear fitting of these data gives 2z/' + 7' = 
98. 1648mc — 37.8980. Recent Monte Carlo renormalization-group analysis of c? = 3 Ising model 
gives a precise estimate of critical point Uc = 0.412051 ± 0.000006 |jl2[, and the series analysis 
of the high-temperature expansions for susceptibility, spontaneous magnetization and specific 
heat gives Uc which is consistent with this |13[. Using this value of the critical point, we can 
read from fig. |l|that 2z/' + 7' = 2.5509 ±0.010, where the error comes from the statistical error 
in the linear fitting. The data given by the approximants with N > 6 also fit to almost the 
same line, but their deviation from the line is larger. 

We plot in fig. |^ the estimate of 2z/' + 7' obtained by fitting the data from all the ap- 
proximants for each fixed L + M + N with N = (p, 0, 1, ■ ■ ■ , 6. We note that the esti- 
mate appears stable for L + M + N > 16. Fitting the data from all the approximants with 
16 < L + M + A^ < 19 and A^ = 0, 0, 1, ■ ■ ■ , 6 gives 



2z/' + 7' = 2.545 ± 0.012. 



(16) 




2.48 
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Figure 2: Estimate of 2v' + 7' by fitting the data from all the unbiased approximants with 
N = (j), 0, 1, ■ ■ ■ , 6 for each fixed L + M + N ; the estimates denoted by (a) and (b) are from 
the high-temperature series and (c) from the e-expansion, respectively. 
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Figure 3: Histogram of the number of biased estimates of 2z/' + 7' for H2; estimates are 



obtained from the approximants with N 



^, 0, 1, 



, 6 and 16 < L + M + AT < 19. 



We also obtain biased estimates by replacing Qm{u) with {u — Uc)Qm{u) in eq. ([151) . The 
estimates of 2z/' + 7' show a good accumulation around 2.55. We show in fig. |] a histogram 
of the number of estimates given by the biased approximants with 16 < L + +A^ < 19 and 
A^ = (j), 0, 1, ■■■, 6. The average and one standard deviation of the data that satisfy 
|2z/' + 7'-2.55| < 0.05 are 

2i/' + 7' = 2.550 ±0.018. (17) 
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Figure 4: Unbiased estimates of 2z/' versus Uc for A2; the data includes the estimates from 
all the approximants with N = (p and 17< M + L < 20. 

Next we give the result of the analysis for A2. The estimates of its critical exponent 2z/' by 
inhomogeneous differential approximants show less convergent results, except for the case of 
N = (f), that is, Dlog Fade approximants. In fig. | we show the plot of the exponent 2z/' against 
the critical point Uc obtained from Dlog Fade approximants of A2 with 17 < L + M < 20. 
Linear fitting of these data gives 2z/' = 43.360'Uc — 16.595. We can read from fig. ^ using 
Uc = 0.412051 that 

2z/' = 1.272 ± 0.004. (18) 

We also list in table |^ biased estimates of the critical exponent 2z/' by Dlog Fade approximants, 
and plot in fig. ^ the estimate of 2z/' by Dlog Fade approximants for each L + M , where we 
have excluded the data that do not satisfy |2z/' — 1.27| < 0.05. All the Dlog Fade approximants 
with 17 < L + M < 20 and |2z/' - 1.27| < 0.05 give 



2z/' = 1.276 ±0.012. 



(19) 



Table 2: Biased estimates of 2z/' from A2; estimates are obtained by Dlog Fade approximants 
with ir = M + L = 17, 18, 19, 20. 



[L/M] 


K=17 


K=1S 


K = 19 


7^ = 20 


[i^-1/1] 


0.4495 


-0.4837 


0.2279 


0.1687 


[K - 2/2] 


-0.4030 


0.8627 


2.1319 


0.8284 


[K - 3/3] 


1.2749 


1.2699 


1.2794 


1.2625 


[K - 4/4] 


1.2751 


1.2737 


1.2731 


1.2579 


[K - 5/5] 


1.2709 


1.2702 


1.2731 


1.2736 


[K - 6/6] 


1.2707 


1.2665 


1.2708 


1.2480 


[K - 7/7] 


1.4254 


1.2698 


1.2707 


1.2716 


[K - 8/8] 


1.2782 


1.2766 


1.2815 


1.2861 


[K - 9/9] 


1.2740 


1.2761 


1.2778 


1.2914 


[K - 10/10] 


1.2655 


1.2769 


1.2735 


1.2807 


[is: -11/11] 


1.2727 


1.2714 


1.2488 


1.2968 


[K - 12/12] 


1.2738 


1.2704 


1.2725 


1.2827 


[K - 13/13] 


1.2688 


1.2722 


1.2735 


1.2745 


[K - 14/14] 


1.2692 


1.2704 


1.2637 


1.2753 


[if - 15/15] 


1.4272 


1.2737 


1.2675 


1.2722 


[A' -16/16] 


1.3291 


1.2950 


1.2889 


1.3058 


[K - 17/17] 




1.2247 


1.2879 


1.2923 


[K - 18/18] 






1.3850 


1.3023 


[ir - 19/19] 








1.2142 



Finally we give the result of the analysis for Ag. The approximants for Ag with A^ = 0, 1,2 
give better converging estimates of 2i/' than the other approximants with A^ > 3, but the 
average of estimates given by these approximants is too small ( 2z/' ^ 1.20 ) and inconsistent 
with the result of the analysis for A2. (If we use the approximants with A^ > 3, we obtain 
smaller estimates of 2i/' ^ 1.17. ) We cannot say that this is due to the fact that the series 
for A2 is shorter than those for A2. The analysis of the series for A2 to order m^^, which is the 
maximum order for Ag, gives more convergent and consistent estimates of 2z/' = 1.28 ± 0.02. 

These values of 2z/' + 7' and 2z/' are to be compared with the result 2z/ + 7 = 2.504(2), 
2.497(4) and 2iy = 1.264(2), 1.260(3) of the high-temperature series [H |ll, 2i/ + 7 = 
2.501(4) and 2z/ = 1.262(3) of the e-expansion p|, and 2z/ = 1.246(6) of the Monte Carlo 



renormalization-group analysis |]I2[. The values of 2z/ + 7 cited here are estimated from the 
values of u and 7 given in the respective references. Our results of 2z/' from A2 [ eqs. (pIS]) and 



(|T9|) ] are not inconsistent with these estimates considering the error bounds. On the other 
hand, our results of 2z/'+7' from /i2 [eqs. (|T6|) and ( [T^ ) ] are 2 per cent larger than the estimates 
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Fig.5 
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Figure 5: Estimate of 2v' by biased Dlog Fade approximants for each fixed L + M ; the esti- 
mates denoted by (a) and (b) are from the high-temperature series, (c) from the e-expansion, 
and (d) from the Monte Carlo renormahzation-group analysis, respectively. 

from the high-temperature series and e-expansion and they are not within error limits. We 
cannot, however, conclude that there is a violation of scaling relation 2z/ + 7 = 2z/' + 7', consid- 
ering the fact that our analysis here does not include the possiblity of confluent singularity. To 
take into account the confluent singularity, we tried an analysis of Roskies-transformed series 



17| with the confluent exponent 0.5, but the result was less convergent. It might suggest 



us to investigate inhomogeneous second-order differential approximants |T8[ , which is another 
method to include the confluent singularity. Longer series might also solve the discrepancy, 
although the estimate from /X2 appears so stable when we change L + M + N from 16 to 19 
in our analysis of inhomogeneous differential approximants as mensioned above ( See fig. 0). 

We find the unphysical but dominant singularity aX u = Ui = —0.2858 ± 0.0003 with a 
critical exponent 

2z/' + 7 (unphysical) = 1.892 ± 0.034, (20) 

from /i2 and aX u = Ui = —0.2858 ± 0.0006 with a critical exponent 



2z/' (unphysical) = 0.811 ± 0.039, 



(21) 



from A2. These values of the critical point are consistent with ui = —0.2853(3) from the 
susceptibility |^, and the critical exponents are so sensitive to the value Mi of the critical 
point if we would try a biased estimation, in which the change of the position of the critical 
point by 0.0005 reduces the critical exponent about 0.08 for fi2 and about 0.03 for A2. 
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5 Summary 



We have calculated low-temperature series for the second moment of the correlation func- 
tion in d = 3 simple-cubic Ising model to order -u^^ by finite-lattice method, from which we 
have obtained the low-temperature series for the second-moment correlation length to order 
M^^ using the known low-temperature series for the susceptibility. The obtained series is 11 
terms longer than those calculated previously. Preliminary analysis of the series by inhomo- 
geneous differential approximants gives critical exponents 2z/' + 7' ^ 2.55 and 2z/' ^ 1.27. 
The latter is not inconsistent with the result from high-temperature series and e-expansion, 
but there is a discrepancy by 2 per cent between the former and the critical exponent from 
high-temperature series and e-expansion. It appears that further analysis of the series should 
be done including the possibility of the confluent singularity. 
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